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■ Abstract 

o 

As a first step to derive the IBM from a microscopic nuclear hamil- 
tonian, we bosonize the pairing hamiltonian in the framework of the 
path integral formalism respecting both the particle number conserva- 
O , tion and the Pauli principle. Special attention is payed to the role of 

the Goldstone bosons. We construct the saddle point expansion which 
reproduces the sector of the spectrum associated to the addition or 
removal of nucleon pairs. 
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c5 ■ 1 Introduction 



The problem of relating the interacting boson model (IBM), so success- 
ful in interpreting the low-energy nuclear phenomenology, to an underlying 
fermionic nuclear hamiltonian has been attacked in a number of papers (see 
in particular [1, 2, 3]). 

In principle, a possible procedure to deal with it could be 

1. to derive an effective interaction in a chosen model space, 

2. to express it in terms of pairing, quadrupole and other forces, 

3. to bosonize the model hamiltonian thus obtained. 
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Less ambitiously one could start by assuming from the outset a multipole 
expansion of the effective interaction in the chosen model space. 

Following these steps one should be able to relate the parameters ap- 
pearing in the IBM to those occurring in the fundamental nucleon-nucleon 
interaction and the bosonic operators to the original fermionic ones. 

A vast literature exists concerning the problem of mapping the fermionic 
model space into the bosonic one [4]. However it appears to us that the 
actual realization of these mappings still deserves futher analyses, especially 
in connection to the nature of the interacting bosons, namely whether they 
are Goldstone bosons or not. Let us comment on why the Goldstone bosons 
should play a role in nuclear structure and how they relate to some of the 
bosons of the IBM. In this connection, we remind that underlying the IBM 
is the recognition that the nuclear interaction is attractive enough to form 
pairs of nucleons. In low energy phenomena these conserve their identity 
inside nuclei, thus becoming the relevant degrees of freedom for the collective 
nuclear levels. The analogy between these pairs of nucleons and the Cooper 
pairs of superconductivity is clearly suggestive and indeed the simplified 
version of the BCS theory of superconductivity represented by the so-called 
pairing hamiltonian was applied to atomic nuclei long ago and was dealt 
with the Bogoliubov-Valatin (BV) transformation. 

The pairing hamiltonian, in its simplest form, governs the dynamics 
of pairs of nucleons moving in a mean field and coupled to zero angular 
momentum (referred to as s-bosons in the IBM). The great success of the 
IBM stemmed from the introduction of pairs of nucleons coupled to angular 
momentum two as well, the d-bosons, which can also be treated through 
a generalized BV transformation. In all cases BV does not conserve the 
number of particles. Although elaborated methods have lately been devised 
to preserve the number of particles [4], it is of importance to realize that this 
feature of the BV transformation is connected to the central concept lying 
at the core of superconductivity, namely the spontaneous breaking of the 
global gauge invar iance related to the particle number conservation. Thus 
in a superconducting system there must exist an associated Goldstone boson. 
Now in the IBM nuclei are indeed viewed as superconducting systems, even 
if only approximately, since they are finite systems. As a consequence, the 
IBM should embody at least one Goldstone boson. Other Goldstone bosons 
should of course appear if other symmetries, like the rotational one, are 
(approximately) spontaneously broken. 

In infinite systems the Goldstone fields carry specific signatures: indeed 
they live in the coset space of the spontaneously broken group ^7(1) of the 
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global gauge transformation e^i^/^^ q being the electric charge, with respect 
to the unbroken subgroup Z2 consisting of the gauge transformations asso- 
ciated with A = and A = nh/q and display only derivative interactions. 
If these distinctive features survive in finite systems, the identification of 
Goldstone bosons becomes of relevance not only for a deeper understand- 
ing of the bosonization mechanism, but also for a convenient choice of the 
variables. One could indeed choose bosonic fields not living in the coset 
space, but at the price of rendering the formalism quite cumbersome, as 
exemplified by the chiral physics [5]. 

Motivated by the above heuristic considerations, in this paper we present 
a new investigation of the program outlined at the beginning of this Intro- 
duction, limiting ourselves to confront the third point, namely to consider 
only the pairing interaction. Yet, since, in a model space, the pairing inter- 
action is an important component of any realistic effective interaction [6], 
then our work may also be viewed as a first step in the derivation of the 
IBM. 

We will use the path integral formalism where the problem of supercon- 
ductivity is readily solved independently from which quantum numbers the 
Cooper pairs carry [5]. The basic question we address is to which extent 
the features of superconductivity in an infinite system, in particular the sig- 
natures of the Goldstone fields, survive in a finite system. Our main result 
is that indeed these signatures are exactly preserved in the pairing model. 
Accordingly we argue that the extension of the latter to include higher mul- 
tipolarities of the force should be easily feasible, as confirmed by recent work 
[7]. 

The plan of the present paper is as follows. 

In Section 2 we discuss from the point of view illustrated above the well 

known spectrum of the pairing hamiltonian, which is characterized by the 
two quantum numbers n (the number of pairs) and s (the number of broken 
pairs) , naturally related to two types of excitations. Those associated with n 
relate to the addition or removal of a correlated pair of nucleons and can be 
viewed as Goldstone bosons stemming from the (almost) broken conservation 
of particle number. Those related to s, namely to the breaking of pairs 
(seniority excitations), we view as corresponding to the Higgs particles. 

In Section 3 we derive the euclidean path integral formulation and intro- 
duce auxiliary bosonic fields via the Hubbard- Stratonovitch transformation. 

In Sections 4 and 5 we set up a saddle point expansion of the effective 
action of the auxiliary fields. The expansion parameter turns out to be the 
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inverse of the energy 

M = f. (1) 

g being the strength of the pairing force and Q = j + 1/2 the pair degen- 
eracy of a single-particle level with angular momentum j (for the sake of 
simplicity we shall consider the case of one level only) . Within such an ex- 
pansion we succeed in reproducing the excitation energies associated with 
the addition or removal of pairs of nucleons. We do not explore however the 
seniority excitations, which are not easily accomodated in the framework of 
our expansion since their energies are of the order of M. A reasonable esti- 
mate, based on the available nuclear phenomenology [8], yields M >3 MeV 
in the region of the Sn isotopes. As anticipated, our analysis shows that the 
pairing model encodes in a striking way the basic features of spontaneous 
symmetry breaking in an infinite system: indeed the field which describes 
the Goldstone excitations lives in the coset of the U(l) symmetry related to 
particle number conservation with respect to the conserved Zi subgroup and 
displays only derivative couplings. 

Then we investigate two paths for selecting a sector of a given number 
of nucleons: one based on the chemical potential (Section 6) and the other 
on a projection operator (Section 7). Since at zero temperature, even in 
the presence of the chemical potential, the number of particles does not 
fluctuate, the two formalisms lead exactly to the same result. 

In Section 8 we derive the hamiltonian of the s-bosons. In the present 
case, since the spectrum of the pairing model is known, this hamiltonian has 
been determined in a simpler, direct way by Talmi [16]. This work, however, 
ignores the Pauli principle, which, therefore, must be added a posteriori as 
an ad hoc prescription. Moreover our approach has the merit that it can be 
applied as well when forces of higher multipolarity are active, thus opening 
the way to the microscopical derivation of the IBM. In Section 9 we present 
our conclusions and outlooks. 

For the sake of completeness, in concluding this Introduction, we mention 
two approaches to the problem of bosonization attempted in the past. 

The first, also limited to the pairing interaction, was based on the use 
of even Grassmann variables in the generating functional [10, 11], but sub- 
stantial difficulties were met in its extension to include the quadrupole in- 
teraction. 

The second was carried out by Mukherjee and Nambu [12] who explored 
in depth the connection between the BCS theory of superconductivity and 
the IBM. These authors, linearizing in the frame of the mean field approach 



4 



a nuclear BCS hamiltonian embodying a contact four-nucleon interaction, 
accounted for second order corrections. They were thus able to derive a 
bosonic Hamiltonian expressed as a sum of Casimir operators, hence qual- 
itatively of the IBM type. However, these authors actually explored an 
infinite homogeneous system and dealt with finite nuclei only in some ap- 
proximate schemes, thus failing to enforce the particle number conservation. 
Moreover their approach appears hardly suitable for a realistic derivation 
of the IBM model, in particular as far as the fermion-boson mapping is 
concerned. 

2 The excitations of the pairing hamiltonian as 
Goldstone and Higgs bosons 

According to the framework outlined in the Introduction we consider the 
schematic pairing Hamiltonian, which accounts for much of the physics of the 
low energy spectra of nuclei. Our aim is to establish a connection between its 
well-known spectrum and the (almost) spontaneous breaking of the particle 
number conservation, which entails the existence of Goldstone and Higgs 
modes. 

In its simplest version the pairing hamiltonian describes a system of 
interacting bounded identical nucleons living in one single particle level of 
angular momentum j and reads 

H= ^2 emXlXm-gAU (2) 

m=-j 

where 

i = ^(-ly— A_^A^, (3) 

m>0 

A^, A are the usual creation and annihilation nucleon operators, m is the 
third component of the angular momentum j, the are the (negative) 
single-particle energies and g is the strength of the pairing force. For the 
sake of simplicity we set 6^ = 6 independent of m. In the Conclusions we 
will mention how and when the level dispersion can be accounted for [13]. 

In this paper we shall consider an even number of identical nucleons only. 
In such a case the energy spectrum is given by the well-known formula [14] 

En,s = 2en - gn{U -n+1) + gs{'D. - s + 1) , n> s (4) 
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where n is the number of pairs and s the pair seniority *. Clearly (4) holds 
valid for 

n<n (5) 

and not only for s < n, but for 

s<n-n (6) 

as well, as a consequence of the Pauli principle. 

In an infinite system the energy of the Goldstone bosons vanishes with 
the associated quantum number. This does not occur in a finite system, but, 
to the extent that the energy spectrum of the latter displays a pattern similar 
to that of an infinite system, it should exhibit two quite different energy 
scales. Actually the excitation energies associated with both the quantum 
numbers n and s appear to be of order gO.. However the Goldstone nature 
of the energy spectrum associated with the quantum number n is clearly 
apparent when one considers the excitations with respect to the minimum. 
This, if n is viewed as a continuous variable, occurs for 

^, = hn + i)-^ (7) 

^ 9 

In fact, since n assumes only integer values, the minimum of (4) takes place 
for 

riQ = [vq], (8) 
[...] meaning integral part. Introducing then the shifted quantum number 

V = n — riQ (9) 

(4) becomes 

Eno+v,s = + "^gvino - vq) - gno{2iyo - no) + gs{n - s + 1) . (10) 

Now we see that the addition (or removal) of one pair of micleons with 
respect to the ground state requires an energy of order g: this is the energy 
of the Goldstone boson. Instead the energy required to break a pair, the 
seniority energy, is of order g fl: this is the energy of a Higgs boson. 

* According to our definition the pair seniority quantum number s is half the usual 

seniority o and corresponds to the number of pairs not coupled to angular momentum 
J = 0, and, as such, blind to the action of the pairing force. 
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We should now point out that using the physical values for fi, g and 
e appropriate, for example, to the Sn nucleus (these can be taken from 
[8], [9]), one would obtain a value for vq corresponding to an unphysical 
nucleus. However the excitation energies of both modes, measured with 
respect to the minimum (see eq. (10)) are essentially independent of the 
single-particle energy e. Hence our argument, although heuristic, remains 
qualitatively correct. The real justification of its validity will be given in 
the next Section. 

Finally, since in the pairing hamiltonian the degeneracy is fixed by 
the model space, the excitation energies related to the quantum numbers v 
and s can be predicted once two conditions are chosen in order to fix the 
parameters e and g. 



3 The generating functional 

As is well-known [15] the path integral must be evaluated in its discretized 
form. The discretized euclidean action of our system is 

j 

-gA{t)Ait - 1) + J2 [Xmit){vt + e)\m{t-l) 
-j 

(11) 

where r is the time spacing, Nq the number of points on the time lattice, 
(v±/)(i) = ±i[/(i±l)-/(t)] (12) 

and 

' -s 




I 



dXdX 



(13) 



the generating functional. Moreover A, A, ^ = J2m>o{~^y~"^^—m^m and 
A = J2m>oi~^y~"^^rnX-m are Grassmann variables. We remind that the 
fermion fields must satisfy antiperiodic boundary conditions in time. 

Now, to cast the action in a form convenient for the saddle point expan- 
sion, we perform a number of manipulations: their role will be illustrated 
when appropriate. First we shift the time label in the variables A (but not 
A) according to 

A(t-l)^A(t) (14) 
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in order to have A and A with the same time argument. This yields for the 
action the expression 

-9A{t)A{t) + E [^m{t) (V+ + e) Xm{t)] \ . (15) 

m=-j ) 

Carrying out next the Hubbard-Stratonovitch transformation, we get the 
new action 

No/2-1 




No/2-1 r 

t=-No/2 I 

+ E [Am(t)(v+ + e)A^(t)] 1. (16) 



m=-j 



Clearly the auxiliary fields rj and rj should satisfy periodic boundary condi- 
tions. 

Finally we introduce the Goldstone field 6 through the polar representa- 
tion [5] for the auxiliary fields 

rj = ^pe^'\ ^ = ^^e-2^^ (17) 

The field p has been placed under square root to avoid the jacobian which 
would otherwise appear. Notice that for this change of variable to be one to 
one (with the only exception of the point p = 0), 9 must vary in the range 
< 6 < TT. Hence the field 9 lives in the coset space of the (almost) broken 
symmetry group U{1) of particle conservation with respect to the unbroken 
subgroup Z2, as appropriate to a Goldstone field [5]. From the periodic 
boundary conditions for the 77-field periodic boundary conditions for p and 
9 follow as well. 

We note now that, after the transformation (17), the 9 field appears in 
the action (16) with non-derivative couplings whereas the Goldstone field 
should display only derivative couplings. However the former can be elimi- 
nated introducing the following transformation on the nucleon fields: 

A™ = e*Vm, A„ = e-*Vm- (18) 
As a consequence of the above transformation the following operators 
g± = ex.p{Ti9)Vf exp(±z^) ± e , (19) 
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whose matrix elements read 
1 



(«")tlt2 



exp < IT 



H2,tl+1 



^ht2 - exp < IT 



t2 



'tlt2 



^t2,tl-l 



+ 6(5, 



tlt2 



eS, 



tlt2 ' 



(20) 

(21) 



will appear in the action. The 6 field appears only in these operators and 
therefore under derivative, as appropriate to a Goldstone field. 
Hence (16) can be recast in the form 



:<in 



No/2-1 r _ _ 

-Na/2 I m>0 



= r 



t=-No/2 
No/2-1 



t=-No/2 



m>Q 



(q^f = -q- 

we find for the fermionic functional integration the result 



Det (- 



q q +9 P] 



(23) 



(24) 



Thus, disregarding here and in the following all the field independent factors, 
we get for the generating functional the following expression 



[dp] / [^6*] exp(-5e//) . 



with 



Seff = r ^ 5p - IV In [-q-q+ + g^p 



(25) 
(26) 



(22) 



X (V'm+5^/p(-l)'■~"('?+)-V-J - g'pi^-miq^r'i^-m + ^-mq^i^- 

Now we first integrate over the fermionic fields and ipm for a given 
positive m: this yields Det(g"'"), independent of m. Likewise, performing 
the integration over ipj^ and V'm for a- given negative m, we get Det(g+ + 
g'^ p[{q^)^^]^), T meaning the transpose operation, again m- independent. 
Lumping the two results together and exploiting the relation 
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Note that the argument of the logarithm is symmetric, but not hermitian. 
The trace must be taken over the quantum number m > and the time. 
The U{1) symmetry is now non hnearly reahzed in the invariance of Sgff 
under the substitution 

e^e + a, (27) 

with a time independent. 



4 The saddle point 



In this Section we look for a minimum of Sgff at constant fields: hence only 
the time-independent component of the p field, to be referred to as p, will 
enter into the effective action. We start by defining 



(28) 



and 



P-I = -V+V,- + e(v+-Vr)+M2 
= -(l-eT)n + M2, 



(29) 



where □ = V^V^- Notice that er cannot be neglected with respect to one. 

Indeed in our calculations we will first perform the limit A^o ~^ oo and then 
we shall let r ^ 0. The effective action at constant fields reads then 



(30) 



The trace is conveniently evaluated in the Fourier representation, yielding 



_ Afo/2-l 

'S'e// = rN^gp - Q ^ In 

no=-Afo/2 



vr 



4(1 - er) sin^ --(no + 1/2) + t^M^ 
No 



(31) 



where the antiperiodic boundary conditions of the nucleon fields have been 
taken into account. Converting the sum into an integral we get 



Seff = iVorj^p-^ln ^ (^rM + ^4(1 -er)+r2M2) | 
= NoT{gp + ne-nM) + 0{T^). 



(32) 
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Notice that the piece fie stems from the term re in (31). The minimum of 
Seff occurs for 

M = f, (33) 

which is independent of e and p. Inserting the above into (28) one gets for 
the value Pq of p at the minimum the expression 

so that Seff at the minimum is 

S, = NoT(^-^-^ + ne^ . (35) 

Although the values for Q, g and e appropriate for the Sn isotopes would 
lead to a negative Pqj actually in selecting a given nucleus (see Section 6) 
the replacement e ^ e = e — n (fi being the chemical potential) should be 
performed. When this is done, as it will be shown in Section 6, pg turns 
out indeed to be positive, as it should be, in the physical range 1 < n < 
attaining its maximum value for n ~ 0,/2, namely for a half-filled level. 



5 The saddle point expansion 

To perform this expansion we start by defining the fluctuation of the p-field 
according to 

p = PQ + r = po(l + ^) (36) 
and by noticing that the generating functional (25) now reads 

/oo rn 
[dr] / [d9]eM-Seff). (37) 
-Po "'O 

Now two cases should be considered: 

a) Pq sufficiently large: then the functional integral defining Z becomes 

gaussian and an expansion in r/pg can clearly be performed. Actually, 
as it will be later shown (sec formula (57) below), pg is indeed large 
when n ~ ^2/2, namely when the level where the pairs live is far from 
being fully occupied or empty; 
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b) po small, which occurs for n ~ 1 or n ~ (see again formula (57) 
below) . In this case the shift in (36) is absent and the p field acts only 
through its fluctuations, which are small, thus assuring the validity of 
the expansion. 

To proceed further we rewrite S^ff in the form 

%/=^E5(Po + + 1^1nP-TVln[l + P(i?i + i?2)] , (38) 
t 

where 

Ri = -q-q" + (V+ + e) (V^" - e) (39) 

and 

R2=g\. (40) 

We set then 

oo 

^e// = E^- (41) 

r=0 

the term Sq being the saddle point contribution, given by (35). This grows 
like Ct^, however it contains also a term of order and a term of order one, 
which should be kept if an expansion in powers of l/fi is sought for. It 
seems to us, however, more convenient to stick to the definition (29) for the 
operator P and to compute the further contributions to the expansion (41) 
(the quantum fluctuations) by developing the logarithm (38): the terms thus 
obtained are naturally organized in powers of M~^. It is worth reminding 
that this expansion does not break the U{1) invariance. 

In the following we shall confine ourselves to evaluate, in addition to the 
first order terms, those quadratic in VtO and r. 



5.1 First order contributions 



These contributions stem from the term linear in r and from the first term 
in the expansion of the logarithm, hence 



-51 =r5^rt-TV[P(Pi+P2)] • 



(42) 



The explicit computation of the second term on the rhs of the above yields 



-Tr (PRi) 



^2 



,2i((?t+i-(?t) 



A9t+i-dt) 



— + 2(l-er) 



(43) 
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where Pu is found to read 



Pft.= 



2M 



1 - er + 



-1/2 



(44) 



By expanding the exponentials up to second order in we get 



-TV(Pi?i) = 



2M 



1 + ( M + -e I T 



■E 



(0. 



t+1 



+ 0{t^) . (45) 



The contribution arising from the third term on the rhs of (42) turns out 
to be 



-lh:{PR2) = -ng^^Pun. 



(46) 



Notably this contribution, linear in r, is canceled by the first term in Si, 
owing to the equation for the minimum of the action S^ff. The cancellation 
holds to the order O(t^), which is the approximation we keep in our analysis 
and in obtaining the equation for the action minimum. 

In conclusion, for the first order contribution to the action we get 



Si = -TriPRi] 



l + (M + -e)T 



(47) 



t=—oo 



where J2t ^"^^s from — oo to oo since we let A?o — oo in evaluating P. We 
note that these contributions are of order J7 and 1, namely are 0(l/il) with 
respect to those of the saddle point. 



5.2 Second order contributions 

We have seen in the previous subsection that all the terms linear in r cancel 
out: hence the r-integration remains undefined. Our aim now is to ascertain 

whether the surviving terms in r stabilize the action. 

Among these we consider the contributions arising from the second term 
in the expansion of the logarithm. They read 

S2 = ^Tr {PRif + TV {PR1PR2) + ^TV (PPa)' • (48) 

In the above the first term is r-independent, the second is linear in r and 
the third one is quadratic. Therefore, for the present purpose, it is sufficient 
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to evaluate the latter. For this we have found 

\T.{PR,f = '^Y.n[Pmfn,. (49) 

Ml 

Hence the integral over r is well defined. 

In conclusion we remark that, as it will be seen in the following Sections, 
in order to obtain the Goldstone bosons energies we must find out how 
they depend upon the single particle energy e. For this purpose we have to 
perform in the integral expressing the generating functional Zi (associated 
with the action ^i) the ^-integration, which appears to be gaussian, but 
actually it is not, because 9 is compact. Yet we can choose Vj0 as a new 
integration variable, thus rendering the integral gaussian. We then get 

1 3 M 3e + gO 

■lnZi = -e + — = . (50) 



A^or 4 2 4 

We notice that this contribution stems from the term (M + 3e/2)r in 
S\, which is irrelevant because it vanishes in the formal continuum limit. 

6 Fixing the particle number by the chemical po- 
tential 

In this Section we apply the saddle point expansion to a specific nucleus 
using the method of the chemical potential. For this purpose we replace e 
with 

e = e — /X , (51) 
/i being the chemical potential. Its value is fixed according to 

where 

N= >^l^rn (53) 

m=-j 

is the particle number operator. Since, however, we shall let A^o ~^ oo 
with r constant, which corresponds to the limit of vanishing temperature 
T = 1/{Nqt), we are allowed to replace < N > with 2n. We also notice 
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that because M does not depend upon e (see eq. (1)), it does not depend 
on /X either. So Eq. (52) becomes 

which gives 

= g{n-n/2-3/8) +e (55) 

for the chemical potential. Hence, in the presence of the chemical potential, 
the energy of the system becomes 

En,o = ^(5o-lnZi) + 2/xn 

= 2en - gn{n -n + 3/4) + f f 5r2 + ^) . (56) 



We thus see from the above that in our approach the excitation spectrum 
of the pairing hamiltonian is reproduced with good accuracy. On the other 
hand the ground state energy differs from the exact value —g[{Q,+l)/2—e/g]^ 
by the quantity {SgCl + g — e)/4, which corresponds to a relative error of 
order 1/0,. 

We conclude this Section by further examining the issue, already ad- 
dressed in the beginning of Section 5, of the validity of our expansion. For 
this purpose it is of importance to assess the size of pg. To this aim we 
replace in (34) e by e and use (55), dropping the term -3/8 in the round 
brackets on the RHS, thus getting 

Pq = n (O — n) . (57) 

Now, when n 0,/2, then the single particle energy e almost coincides with 
the chemical potential p. In such a situation e is almost vanishing and, from 
(57), po — ^^^/4. This large value corresponds to the situation when the 
level where the pairs live is neither fully filled nor almost empty. On the 
other hand pg attains its lowest value when n = 1 or n = fi. It is remarkable 
that even in these cases, where an expansion in 1/po cannot be performed, 
our approach still yields the correct excitation spectrum of the system. 



7 Fixing the particle number by the projection op- 
erator 

Owing to the importance of properly fixing the particle number n, in this 
Section we address the problem through an alternative procedure, namely 
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by introducing in the path integral the projection operator 



(58) 



Using then the variables (18) and performing the Hubbard-Stratonovitch 
transformation as previously done, we get for the generating functional the 
expression 



da 
-TT 27r 



-5(") 



(59) 



where 



No/2-1 ( _ _ 

t=-No/2 I m>0 
+9VP {^m'>P-m + i^-mipm) \ - 2NoT{a - e)n 



(60) 



(remember that the label n indicates the number of pairs), having defined 



a = e + 



Not 



and 



qf = exp{^ie)Vt exp{±i9) ± a . 



(61) 



(62) 



Next we carry out the integration over the fermionic degrees of freedom, 
getting for the partition function the expression 



cx j 2^ da J [dpde]exp[-S^J}\{a,p,9)] , 



being 



qin) 
^eff 



-2NoT{a - e)n + Tj29P-^^^ -Qalt+d^P 
t 

At constant fields (64) simplifies to 

gp + 2ne + (O - 2n)a - n^a'^ + g'^p 



S^''\a,p) = Not 



(63) 



(64) 



(65) 
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which is stationary when 



^^"'<'-^)-~„.,fi-^^Uo (00) 



9p " ^ V 2v/^^TF^ 

and 

A^oT(^-^ + n-2n) =0. (67) 



da 

The solutions of the above equations read 



Po = n{n-n) (68) 

and 

(70 = |(J^-2n). (69) 
Finally the effective action (65) at the minimum ctq turns out to be 

S^"^ ((70, po) = NqtE^''^ = NoT[2ne - gn{n - n)] , (70) 

which differs from the zero seniority spectrum of the pairing hamiltonian. 
Indeed the latter has + rather than CI, inside the round bracket on the 
RHS of (70). 

Note that the action (65) is an analytic function of a inside an integration 
path deformed to encompass the saddle point o-q. This path goes from 
e — ztt/Nqt to c7o — in/NoT along a straight line parallel to the real axis, 
then from ctq — itt/Nqt to (Tq + itt/Nqt along a straight line parallel to the 
imaginary axis and finally it goes back from ao + iir/NoT to e + itt/Nqt. 
When iVo — GO, the contributions coming from the paths parallel to the 
real axis cancel each other, while the one parallel to the imaginary axis 
vanishes: hence there are no corrections to the saddle point contribution. 

Next, with the aim of checking the results obtained in the framework of 
the chemical potential method, we evaluate, using the projection operator, 
the first order correction in the saddle point expansion in p and 6. For this 
scope we set, as in (36), 

p = Po + r. (71) 
The action (64) can then be recast as follows 

siffC.S^^\ao,Po) + sl-^ (72) 
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where 

Si'^) = N^rgr - TV[P(ao, Po)(4"^ + 4"^)] (73) 

with 

4"^ - ^o)][qt - (a - ao)] + (V+ + ao)(V,- - ctq) (74) 

and 

4"^= A- (75) 

The contribution 

Tr[P(ao,Po)4"^] =5'TrP(ao,Po)r = /-$]r (76) 

^ t 

cancels the first term in (73). Hence the latter simply becomes 

5f) = -TV[P(ao,Po)4"^]- (77) 
A calculation, similar to the one carried out in Section 5.1, yields then 

Using now the expression (69) for ctq we get 

TY,e{-n)e (79) 
t 

and by performing the ^-integration (again using (69)) we finally obtain the 
first order energy 

^(«) ^ ^(") =2ne- gn(^n-n + ^ + ^gn, (80) 

which coincides with the excitation spectrum obtained with the chemical 
potential. 

8 The hamiltonian of the s-bosons 

To complete our program (restricted, we remind, to the pairing potential) we 
derive below the bosonic hamiltonian corresponding to our effective action. 



-Tr 



PR); 



(n) 



si") = 



1 . /5^ 3 , 
- l + g(-n--n\r 
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The most general, particle conserving, quartic hamiltonian for a system of 
s-bosons, confined lo live in one single particle level, reads in normal form 

H{b\h) = hh^b + vhWhh, (81) 

0^ ,h being bosonic creation-annihilation operators acting in a Fock space. 
They satisfy canonical commutation relations. The values of the parameters 
/i, V were obtained in ref. [16] in such a way to yield the pairing Hamilto- 
nian spectrum. Of course the Hamiltonian thus obtained, being intrinsically 
bosonic, patently violates the Pauli principle and therefore the condition 
n <Q. should be added a posteriori^ when using (81) in describing a system 
of fermions. Here, to show how this condition naturally emerges instead in 
our framework, we obtain the parameters h and v with our methods. To 
this purpose we write the path integral associated to (81), namely 

Z = J[db*db]exp{-S), (82) 

where 

No 

S = tJ2 {K+i'^tbt + H{bl^M] (83) 

and the b* , h are holomorphic variables satisfying periodic boundary condi- 
tions in time. We now compare the above to our effective action. 
To this end we introduce the polar representation 

b = Vpexp(ze), h* = Vpexp(-i0) (84) 

in terms of which the generating functional and the action read 

/•oo 

Z = / [dp] / [d^]exp(-5) (85) 

S = l^pt+i exp(-i0t+i)Vt [v^exp(z6't)] 

t 

+ -H' [VPt+iexp(-i6't+i),^exp(i6't] |. (86) 

Now we again look for a minimum of S at constant fields and perform a 
saddle point expansion. The calculation is basically the same as the one 
previously developed and hence will not be reported here. We only remind 
that such an expansion holds valid only for n < $7, as it follows from the 
positivity of po, given in eq. (57): thus in our approach such a condition, far 
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from being artificial, is necessarily implied by the formalism itself. We only 
quote the result of the comparison with our effective action: it yields 



Obviously the considerations following (41) hold valid as well here. Thus h 
will be affected by an error of order 1/M. Indeed for the Hamiltonian (81) 
to reproduce exactly the pairing spectrum it must be 



Finally, it is of importance to stress once more that the above discus- 
sion can be generalized to include other types of bosons, as, for example, 
those appearing in the quartic hamiltonian of the IBM model, respecting 
basic symmetries like particle number conservation, rotational invariance, 
etc. It is then clear that our approach opens the way to microscopically 
deduce the Arima-Iachello model: in this case to fix the coupling constants 
one should again write the corresponding path integral and compare the re- 
sulting bosonic action with the effective one found extending the procedure 
developed in this paper to a fermionic hamiltonian including forces of higher 
multipolarity. 

9 Conclusions 

In this paper we have carried out an investigation concerning the possibil- 
ity of a systematic bosonization of a realistic nuclear hamiltonian for the 
description of the low-lying sector of the nuclear spectrum. Our study is 
admittedly preliminary since it is limited to the pairing interaction, which 
is of course only a component (although an important one) of an effective 
interaction pretending to be realistic. However, having overcome the main 
difficulty we expected, namely the one of going from an infinite to a finite 
system, we are now confident to be able to solve the bosonization problem 
in the presence of other types of interaction [7] . 

Our approach is based on the concept of symmetry breaking and on the 
related properties of the interaction among the bosonic fields: indeed this 
framework is the most suitable for deriving, rather than assuming, a model 
like the IBM one. 

To develop our scheme we have used the path integral formalism because 
of the large flexibility it allows both in choosing and in dealing with the vari- 
ables appropriate to the problem. This has lead us to deduce an asymptotic 
expansion for the system's spectrum in the parameter = 



h = 2e- gQ- g/A, v = g. 



(87) 



h = 2e — gQ, v = g. 



(88) 
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The euclidean path integral clearly encompasses the whole Fock space 
of the system. To deal with a specific nucleus a given number of pairs must 
be selected (a procedure not to be confused, of course, with the projection 
of the particle number when this is violated). In our approach a definite 
particle sector has been chosen using both the chemical potential and the 
projection operator methods: the two turn out to be completely equivalent. 

Notably our expansion to first order reproduces with good accuracy the 
energy of the pair addition and removal modes (or, in the language of the 
IBM, of the s-bosons). Moreover the requirement that should be pos- 
itive entails the inequality n < Q (see again eq. (57)), thus implementing 
the action of the Pauli principle. This is conceptually important for the 
consistency of our scheme, which respects both the Pauli principle and the 
particle number conservation. Worth noticing is that this crucial feature is 
absent in the framework developed in [16, 17], but correctly dealt with in 
many treatments based on the mapping procedure [1, 4]. 

Our expansion cannot account for the seniority excitations, whose ener- 
gies, in the framework of the pairing Hamiltonian (2), are larger than M: 
hence they should be separately treated. For this scope clearly an analysis 
of the excitations associated with the field p should be performed. In this 
connection it is worth reminding that in the scheme of the BV transforma- 
tion these modes are described in terms of quasi-particles whose energies, 
as well-known, are expressed in terms of the gap A. In an infinite system 
A is associated with the order parameter by setting a non-zero vacuum ex- 
pectation value for the pair (Cooper) field and it signals the onset of the 
superconducting phase; in a finite system it measures the strength of the 
mean field (referred to as the "pair field") which, in the BV scheme, lin- 
earizes the pair interaction. 

Finally we should observe that in our analysis the impact on the exci- 
tation spectrum stemming from the removal of the degeneracy of the levels 
where the pairs are sitting and of the higher order terms in the saddle point 
expansion has not been explored. Concerning the first issue, if the spacing 
between the single particle energy levels is small with respect to M, it can 
easily be accounted for within the present perturbative scheme. Indeed this 
approach has already been pursued in ref. [13]. Actually the Sn isotopes are 
well suited for such perturbative treatment since here the distance between 
the single particle energies of concern appears to be much smaller than M. 

Concerning the second point, we can only say that it would certainly be 
both interesting and important to examine in more depth our asymptotic 
expansion. 
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